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Abstract 

In this note, we prove that the space of all admissible piecewise linear metrics 
parameterized by the square of length on a triangulated manifold is a convex cone. We 
further study Regge’s Einstein-Hilbert action and give a more reasonable definition of 
discrete Einstein metric than the former version in Finally, we introduce a discrete 
Ricci flow for three dimensional triangulated manifolds, which is closely related to the 
existence of discrete Einstein metrics. 


1 The space of piecewise linear metrics 

Consider an n dimensional compact manifold M with a triangulation T. The triangulation 
is written as T = {7o,71,-- - ,7^}, where 7i (0 < i < n) represents the set of all i 
dimensional simplices. A piecewise linear metric is a map Z : 7i ^ (0, +oo) making each 
simplex an Euclidean simplex. 

There are two disadvantages to think of I as the analogue of smooth Riemannian metric 
tensor g. For one thing, we know that OJIt) the space of all admissible piecewise linear 
metrics, is not convex (although it is a simply connected open set). For another, the 
scaling property of I is not good enough. If the smooth Riemannian metric tensor g scales 
to eg in the smooth manifold M”, then the length l{'y) of a curve 7 : [0,1] ^ M scales to 

VcKl)- 

If we take as the direct analogue of metric tensor g, both the above two disadvantages 
can be overcome. The idea of considering the square of Z, not Z itself, as an analogue of 
smooth Riemannian metric tensor comes naturally from the former work by the first 
author and Xu [1] , where the idea has been used for piecewise linear manifolds with circle 
or sphere packing metrics. Firstly, we have 

Theorem 1.1. For manifold M" with triangulation T, denote gij = ifj for each adjacent 
edge i ^ j. Then the space of all admissible piecewise linear metrics parameterized 
by gij, is a nonempty connected open convex cone. 
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Proof. Rivin im first observed this fact for a single simplex case. Gu et al [8] proved 
this fact for re = 2 by direct calculation. The proof here follows from Rivin’s idea. For an 
re-simplex A embedded in the Euclidean space, we label all vertices as vo,vi, - ■ ■ ,Vn and 
all edges as Iqi, ■ ■ ■ ,ln-in- For brevity, let re* = idAhl)^ then we need to show 

= ,iL G I ^01) • ■ • Jn-in are edges of some Euclidean re-simplex} 


is convex. Construct a map from to the set of all symmetric re x re matrices, which 
transforms • • • , to 
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The above matrix is the Gram matrix of re linear independent vectors 01,02, • • • ,0re and 
hence is positive definite. Obviously, the map is injective and surjective. Note that the 
set of all positive definite re x re matrices is a nonempty open convex subset of R"' . Thus 
is also a nonempty open convex subset of R” . 

Next we prove the theorem for general triangulations. Assuming all edges are labeled 
in turn as ei,--- ,6^, where m = |Tl|. Then for any re-simplex A = (uq,--- ,Un) with 
edges eii,...,ei^*, (ii,--- ,in* G {1,2,-•• ,m}), denote 




we have 


ml= n ml. 

This implies that 501^ is a nonempty connected open convex cone of R™. 


□ 


2 An interpretation of Regge’s Einstein-Hilbert action 

Now we deal with three dimensional case. We give an interpretation to three dimensional 
Regge’s Einstein-Hilbert action. The idea here is nature when taking 1“^ as the analog of 
Riemann metric tensor g. Given an Euclidean tetrahedron {z,j. A:,/} G Ta, the dihedral 
angle at edge {i,j} is denoted as The discrete Ricci curvature Rij at edge {i,j} is 

Rij = 27r - ^ (2.1) 

{i,j,k,l}£T3 
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where the sum is taken over all tetrahedrons having {i,j} as one of their edges (If this 
edge is at the boundary of the triangulation, then the discrete Ricci curvature should be 
Rij = TT — j f. z }675 Pij,ki)- Regge’s Einstein-Hilbert action [TOl [U [7] in this case is 

£ = ^R^jhj. ( 2 . 2 ) 

Denote ai^jki be the solid angle at vertex i in a single tetrahedron. Cooper and Rivin 
[2] once defined a combinatorial scalar curvature at vertex i as 

= dvr - ^ 

where the sum is taken over all tetrahedrons having i as one of their vertices. This 
dehnition satishes the following combinatorial equality at each vertex i 

Sf^ = Y,R^R ( 2 . 3 ) 

j-^i 

which can be proved by Euler’s characteristic formula for spheres. However, Cooper and 
Rivin’s curvature is scaling invariant, hence it does not perform as well as smooth scalar 
curvature when the metric scales. In the following, we denote as the inverse of gij, i.e. 
g^^ = g~-^. For any function / defined on all edges, we can take a combinatorial trace at 
vertex i with respect to g as {trgf)i = J2jr^i 9^^ fij- Thus we may define the discrete scalar 
curvature as 5 * = tVgR, i.e., the discrete scalar curvature at vertex i is 

Si = {trgR)i = Y,g'^Rij- (2.4) 

jr^i 

We also define a combinatorial volume element at vertex i as 

= (2.5) 

The volume element and scalar curvature can be multiplied together by a combinatorial 
multiplication *, which is defined as 

{S *V)i = Y, • ( 2 - 6 ) 

jr^i jr^i 

Taking (5 * V)i as a combinatorial analogue of smooth Sdvol, we can get 

£ = \Y.^S*V), = Y,R^3h, ( 2 . 7 ) 

i j^i 

which is a combinatorial analogue of smooth Einstein-Hilbert functional Sdvol. 
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3 Discrete Einstein metric 


Let V = ^ be the discrete volume of the manifold M. It is a combinatorial analogue 

of smooth volume Jj^dvol. The normalized Einstein-Hilbert functional 

fj^Sdvol 

plays an essential role in studying three dimensional Yamabe problem. This inspires us to 
consider a normalized Regge’s Einstein-Hilbert action 


Qia) = 

E3 


i^j 



(3.1) 


Since Rij and gij = tfj are somewhat like smooth Ricci curvature Ric and smooth 
metric g respectively, the Einstein metric g with Ric = Xg on smooth manifolds M can 
be transformed to a discrete version, i.e., a combinatorial metric g with R = on 
a triangulated manifold This fact inspires us to define the following discrete 

Einstein metric 


Definition 3.1. Given a compact manifold with triangulation T. A piecewise linear 
metric g is called a discrete Einstein metric if there exists A such that R = Xg. 


It is easy to see that, if R = Xg, or say s' is a discrete Einstein metric, then 



i^j 


Theorem 3.2. Given a manifold with triangulation T, then a piecewise linear metric 
<7 is a discrete Einstein metric if and only if it is a critical point of the normalized Regge’s 
Einstein-Hilbert action. 


Proof. The Schlafli formula [9] says 

^ ^ ^ijdRij — 0 . 


Hence we get d£ 


Then 


V~3 _1 

^gQ = —^diag 2 {c/i, • • • ,5m}(R - A^), 


which implies the conclusion. 


□ 
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We want to know how many discrete Einstein metrics there are for a fixed three 
manifold M with triangnlation. How to find them? Can we triangnlate M snitably, so 
as M admits a discrete Einstein metric? We introduce a new topological-combinatorial 
invariant, which carries the information of the triangnlation, 

Ym,t = inf Q{g) (3.3) 

associated with a fixed triangnlation T on a fixed manifold . We also introduce a 
pure topological invariant Ym = supg-Y m,t associated with a fixed manifolds M^, where 
the supremum is taken from all triangulations of M. We hope the study of these two 
combinatorial and topological invariants will deepen the understanding of discrete Einstein 
metrics. 


4 Discrete Einstein metric of ct-order 


Similar to the works by Ge and Xu lais], almost all procedures about discrete metric gi 
and discrete curvature Ri can be generalized to a order. We can define a-order discrete 
Einstein metric, i.e. a metric I satisfying R = Aq,/" (a € R). 

The logic behind this procedure is to take as a metric (of a order). From the 
viewpoint of Riemannian geometry, a piecewise flat metric is a singular Riemannian metric 
on M^, which produces conical singularities at all vertices. For any a E R, a metric g 
with conical singularity at a point can be expressed as 


ff(z) 


= 


dzdz 


locally. If letting f{z) = — In(a^), then g{z) = \dz°‘\‘^. Compared with |dz“|, the 
a-metric /“ could be considered as a discrete analogue of conical metric to some extent. 
Inspired by the essential role the conical metric plays in the study of smooth canonical 
or extremal metrics, we want to know what happens if we triangulates the manifold and 
evaluates the piecewise linear metrics of it. It seems that canonical or extremal metrics 
in the piecewise linear metric case should be a metric parallels to the discrete Ricci 
curvature R, i.e., a discrete Einstein metric of a-order. 

It is easy to see that, if Hs a discrete Einstein metric of a-order, then 


^ijhj 

E ' 


(4.1) 
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Furthermore, a-order discrete Einstein metrics are critical points of the a-order normalized 
Regge’s Einstein-Hilbert action 


Qa — 


ir^j 


(4.2) 




By this definition, the discrete Einstein metric in Definition 13. II is in fact a 2-order discrete 
Einstein metric, the discrete Einstein metric defined in [3] is actually a 1-order discrete 
Einstein metric, and the discrete Ricci flat metric is a 0-order discrete Einstein metric. 


5 Discrete Ricci flow 

Definition 5.1. Given a manifold with triangulation T, the discrete Ricci flow is 


dg-. 




= -2Rij, or 


dg 


= -2R. 


dt u; “ 

It is useful to consider the normalized discrete Ricci flow 

dgij 2 dg 2 

— = -rg,,-2R,,, or - = -rg-2R, 


(5.1) 


(5.2) 


where r = 3A = 3£jV. 

Discrete Ricci flow (j5.2|] takes the same form as the smooth Ricci flow ^ = ^rg — 2Ric 
with n = 3. It is easy to have 

Proposition 5.2. Along the normalized discrete Ricci flow (15.211 . R is a constant, £ is 
descending and bounded. Moreover, Ricci flow (15.2p is a negative gradient flow. □ 

Theorem 5.3. If the solution g{t) of the discrete Ricci flow (|5.2p converges to a non¬ 
degenerate metric ^(-l-oo), then g{+oo) is a discrete Einstein metric. □ 


Recall that the Schlafli formula [9j implies ^ hdRi = 0 and hence d£ = Yh Ridl 

i=l i=l 

Thus we have V/T = = R- Denote L = ; then L = Hessi£ is 

symmetric matrix. However, generally L is not positive dehnite [3]. Denote 


Al = 


d{Ri,--- ,Rm) 

d{gi,--- .gm) ’ 


then Al could be thought of as a discrete version of Lichnerowicz Laplacian. Denote 
XinfiA]^) as the smallest eigenvalue of then we have 
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Theorem 5.4. Let g* be a discrete Einstein metric. If Ain/(A2) > A* and the initial 
metric g{0) deviates from g* not so much, then the solution to discrete Ricci flow (|5.2I1 
exists for all f > 0 and converges exponentially fast to the discrete Einstein metric g*. 

Proof. The proof of this theorem is similar to those appeared in BUM- We 
show that g* is an asymptotically stable point of the normalized discrete Ricci flow (15.21) . 
Denote T{g) = ^rg — 2R = 2{Xg — R), then g* is a critical point of T{g). We denote 
S = diag{li, ■ ■ ■ ,lrn} and differentiate r(g() at g*, 

^Dg.T{g) = A*(/™ - - Al. (5.3) 

If Xinfi^l) > A*, then by similar methods in BUM, we can show that all eigenvalues of 
the matrix Dg*T{g) are negative (up to scaling of metric l^). The proof is completed. □ 

Similarly, for any a G IR — {0, —1}, we can also define the normalized a-order discrete 
Ricci flow as 

dl°‘. 

= 2XJfg - 2Rij, (5.4) 

which can be written in vector form as 

dl° 

V = 2A„r - 2R. (5.5) 

dt ^ ^ 

Remark 1. We have recently learned that R. Schrader got the same result (Theorem 3.1, 
m) as our Theorem o in this paper. R. Schrader m also provided some interesting 
analogues of Ricci tensor and Ricci flow in the theory of piecewise linear spaces. 


6 An example: 16-cell triangulation 

Discrete Ricci flow provides an efficient way to find discrete Einstein metrics on trian¬ 
gulated manifolds. To see the benefit of this method, we take the 16-cell triangulation 
of as an example. Here 16-cell is defined as a triangulation 71 of like this: take 
Hi = (1,0,0,0), A2 = (-1,0,0,0), Bi = (0,1,0,0), B 2 = (0,-1,0,0), Ci = (0,0,1,0), 
C 2 = (0,0, —1,0), Di = (0, 0,0,1), D 2 = (0,0,0, —1) as the vertices of 7^; PiQj{{P, Q} G 
{A,B,C,D},i,j = 1,2) as the edges of 7^; PiQjRk {{P,Q,R} C {H, R, C, 71}, i, j. A: = 
1, 2) as the faces of 71; and regular tetrahedrons AiBjCkDi{i, j, k,l = 1, 2) as the tetrahe¬ 
drons of 71. 

Now we consider a topological triangulation T of S^, which has the same combinatorial 
structure with 16-cell 71. Obviously, 7” carries a trivial a-order discrete Einstein metric 
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for each a, i.e. a metric with the same value on all edges. We want to know whether 
there exist other non-trivial discrete Einstein metrics? Let Ip^Qj be the lengths of PiQj, 
and denote gp^Qj = Consider the normalized discrete Ricci flow 


dgPjQj 

dt 


= ‘^^gPiQj - ‘^RPiQj- 


( 6 . 1 ) 


We now specialize the parameters in equation (|6.ip : 

1. A = the sum is taken over all the twenty four edges in T. 

^{PiQj}'-PiQj 

2. To illustrate the calculation of RpiQ^, we take RaiBi as an example. Note that 

RAiBi =277 - I3aiBi,CiDi - /3aiBi,CiD2 - I3AiBi,C2Di - /3aiBi,C2D2, (6.2) 


where l3AiBi,CiDj ihj = 1,2) are the dihedral angles at edge AiBi, 


sin^AiBi.CiDj 


^J-AiBiVAiBiCiDj 
2SAiBiCiSAiBiDj ’ 


(6.3) 


12^2 


^AiBiCiDj —^AiBjCiDji^AiCi + + ^BiCi + ^BiDj ~ ^AiBi ~ ^CiDj) 


^ ^ o J ^ U ^ J ^ ^ ^ J ...... J 

+ ^AiDd'siCii^AiBi + ^CiDj + ^AiCi + ^^BiDj “ ^‘aiDj “ ^IbiCi 
+ ^AiCi^BiDj i^AiBi + %iDj + ^MDj + “ ^Ai 

72 ]2 ]2 i2 i2 j2 

- '■AiBPBiCdAiCi - ‘AiCdAiD/Ci 

_ /2 j‘2 )2 




BiDj. 


j2 j2 j2 

'iDi ~ '■AiBi'-AiDj‘‘BiDj 


(6.4) 


SAiBrCi = \JPAiBiCiiPAiBiCi “ UiBi)(RA iBiQ “ lAiCi){PAiBiCi “ (6.5) 

SaiBiDj = Y^RAiBiDj (RAiBiDj — lAiBi){PAiBiDj — lAiDj){PAiBiDj “ IriDj), 

( 6 . 6 ) 

where PAiBiCi = (^AiSi +UiCi + ^SiCi)/2 and PaiBiDj = (^AiBi +lAiDj +^BiDj)/ 2. 
Similarly, we can calculate all the other Rp^q^. 


It is easy to see that gPiQj{t) = 1 is a trivial solution of equation (j6.1|) . With the aid of 
Matlab software, we can solve equation (16.111 numerically. The numerical calculations show 
that discrete Ricci flow equation (16.ip may develop singularities for some initial values. 
However, if we choose the initial values appropriately, the solution g{t) to equation (16.11) 
can converge and hence deforms the metric g{t) to some discrete Einstein metric g*. 
Write the solution g{t) as a 1 x 24 row vector with order IaiBi^ ^AiB 2 ’ ^A 2 Bi’ ^A 2 B 2 ’ 
^AiCi’-- - >^AiDi>--- JbiCiP-- ^^BiDp--- Denote 1„ as a 1 x n row vector 

with all entries being 1, then we can give two specific discrete Einstein metrics as follows: 







1. g* = (ali2,61i2), where a = 10.0095, b = 6.9633. 

2. g* = (els, dl4, el4, els), where c = 58.7223, d = 64.9735, e = 52.1413. 
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